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Chapter 1 


Introduction 


The objective of determining the constitutive parameters of a material 
placed in a fixture can be accomplished by S parameter measurements 
conducted by broadband network analyzers [1]. A common piece of ap- 
paratus used in the measurements is a waveguide fixture. Traditionally 
when a waveguide is used, its cross section is completely filled with a fixed 
length of the material to be measured. A method of parameter determi- 
nation is now desired where the technique used yields results which is not 
significantly affected when the fixture is subjected to a high temperature 
environment. Such an environment is likely to introduce error to a com- 
pletely filled wave-guide due to temperature expansion of the material and 
the waveguide itself causing air gaps between the fixture and material. Ac- 
curacy may even be diminished when undesired air gaps already exist or 
when there is not a sufficiently tight fit between the waveguide walls and the 
material. Error correction formulas compensating for the air gaps have aj> 
peared previously [2,3,4]. However, one is valid only for low loss dielectrics 
and the other is merely an empirical formula. 

An alternate method is investigated for the rectangular waveguide where 
a known air gap exists between the waveguide walls and material sam- 


1 


pie [5,6,7]. The goal is to minimize the errors which occur due to limita- 
tions in the accuracy of the dimensions of the materia] and change in the 
dimensions caused by thermal expansion. Two different geometries are con- 
sidered. The first geometry has the air gap oriented parallel to the shorter 
wall. The other has the gap running parallel to the longer wall. The main 
disadvantage of these alternate geometries is the increased complexity they 
add to their S parameter expressions, thereby increasing the amount of time 
in extraction of the proper set of constitutive parameters. 

This report will first discuss the traditional method of extracting e and 
/ x from S parameter measurements using fixtures with their sample region 
completely filled [8]. This direct method does not work when the cross- 
sectional region is not completely filled due to the presence of higher order 
modes. An alternative approach is presented to accomplish the e and \i 
parameter extraction process. Unlike the former approach, explicit rela- 
tionships between the material parameters and the S parameters are not 
possible. The latter approach is dependent upon an iterative technique 
(Newton-Raphson) to acquire the material parameters from the S parame- 
ters through the use of an analytical model. 

The presence of partially-filled cross-section requires a detailed field 
analysis. Such a field analysis is described first for two rectangular waveg- 
uide samples geometries followed by a parameter study for the S parameters. 
Next, the actual iteration technique is discussed to obtain the constitutive 
parameters from the measured S parameter data. This technique is applied 
and demonstrated for two different measurement fixture configurations, i.e., 
a one port and two port fixture. The measured S parameters for the differ- 
ent fixture configurations are examined and compared to expected values 
based upon independent extracted values for t and ji. This is followed by 
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a display of extracted parameter values using the developed techniques. 
Finally, the conclusions of these methods for parameter determination and 
suggestions for further exploration are given. 
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Chapter 2 

Material Parameter Extraction 
Using Waveguides 

A common approach for parameter determination is to obtain reflection and 
transmission measurements for a fixed length of material in a waveguide 
fixture (coaxial or rectangular) with the fixture cross section completely 
filled. The parameters extracted from such an approach are subject to error 
from sample condition (homogeneity, dimensions and shape distortions) and 
measurement hardware calibration. Examples of error magnitude will be 
presented for some of the errors. The error which occurs from a partially- 
filled fixture cross section is presented in the following chapter. The basic 
equations used for such parameter extraction are presented here. 

The S parameters Sn and S 2 i can be measured by automated broad- 
band network analyzers in the frequency range where only the dominant 
mode is able to propagate in the unfilled sections of the waveguide. These 
measurements can be related to the reflection (-Sn) and transmission (P) 
coefficients of a material junction of semi-infinite length. The relationships 
are 


(1 - P 2 (£ ) /X))5 1 i(g,/X) 

1 - P 2 (e,n)s 2 u (e,n) 


( 2 . 1 ) 


and 


S 21 — 


(1 - 3 2 n (e,fi))P(e,fi) 

1 - P 2 (e,n)s 2 n {e,fi) ‘ 


The previous expressions can be rewritten to express Sn and P in terms 
of Si 1 and S 21 . These expressions are 

p = . Sl ‘ ~ % r ’f (2.3) 

1 ~ (Sn + SiiJ'Sn 


s n = K± VK 2 - 1 


where 


S 2 — S 2 +1 

zr _ J 11 °21 ' 1 


The sign in Equation (2.4) is chosen to insure |sn| < 1. 

These new characteristics, and P, can be related to the material 
parameters. The relationships for a coaxial waveguide fixture are 

p _ e -i*V^d (2.6) 


J*-- 1 

vF+T' (2 ' 7) 

Note that e = e r e a and /x = fi r fi 0 where the subscript notation refers to 
relative and free space values. 

The corresponding relationships for a rectangular waveguide fixture are 
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p — ly/'rVrd 


( 2 . 8 ) 


and 

( 2 . 9 ) 

where X Q — yj\ — (|^) 2 and Xj = <Jl - (— ^l _ ) 2 with A 0 being the free 
space wavelength, c the speed of light and a is the width of the rectangular 
guide. 

The final expressions for the constitutive parameters using a coaxial 
waveguide are given by 

fir = \/xy ( 2 . 10 ) 



and 



where 



and 



( 2 . 11 ) 


( 2 . 12 ) 


( 2 . 13 ) 


The final expressions for the constitutive parameters using a rectangular 
waveguide are given by 
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(2.14) 


1 + ^11 
1 — 5n AX Q 


and 


_ (± , _J_\ K 

* r \A* + (2a) 2 / fi r 

where 


(2.16) 

Note that In p is modulo j2nn where n is the integer of The proper 
value of n can be estimated with the integer value of 

(2.17) 

where is the derivative (slope) of the phase for the transmission coefficient 
P with respect to frequency and / is frequency. 

An important aspect in parameter determination is the error occur- 
ring from calibration and sample dimensions. The following examples are 
for coaxial and rectangular waveguide fixtures with sample parameters of 
e r = (4.0,— j. 20) and n r = (1.2, — j.12) and a length of d=. 3". Figures 2.1 
through 2.4 illustrate error from the phase reference plane position for phase 
offsets in Sn corresponding to distances of .001” and .003”. The error is 
defined as 


n = Int 


fV 


2n 


1 

A 2 


-L.nl 

2nd P 


i2 


(2.15) 


„ Extracted Value - Exact Value nntw 

% Error = - x 100%. 

Exact Value 


(2.18) 


The phase offset error (position) results in a comparable amount of error for 
both fixture types with the imaginary component of t being very sensitive 
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to the position of the reference plane. Similar results also occur for the 
parameter. 

Figures 2.5 through 2.8 illustrate the error when the sample length is 
assumed to be .3" but is actually .302" for the solid curves and .298" for the 
dashed curves. The fixtures perform differently in this case with the coax- 
ial fixture yielding error which remains essentially constant with varying 
frequency. The rectangular waveguide fixture has error which decreases as 
frequency is increased. For both fixtures the error in extracting constitutive 
parameters is not very sensitive to inaccurate sample lengths. 

Other sources of error arise from other discontinuities occurring in the 
measurement system. The best solution so far to eliminate or minimize 
these errors lies in a calibration technique [9]. It is also advisable to use 
a fixture with sufficient length to separate the returns due to the fixture 
junctions from the sample return in the interior region of the fixture and 
to insure that only the dominant mode is incident on the sample. 
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Figure 2.1: Normalized error for e when S u has a position error of .001" 
(solid) and .003" (dashed). The sample parameters are e T = (4.0, — j.20) 
and fi r = (1.2,— j. 12) with d = .3". Coaxial airline fixture. 
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Figure 2.2: Normalized error for ft when Sn has a p osi tion error of .001" 
(solid) and .003" (dashed). The sample parameters are e r — (4.0, —.7.20) 
and n T — (1.2, —j. 12) with d — .3". Coaxial airline fixture. 
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Figure 2.3: Normalized error for e when S n has a position error of .001" 
(solid) and .003" (dashed). The sample parameters are e r = (4.0, —j. 20) 
and [i r = (1.2, — j.12) with d = .3". Rectangular waveguide fixture. 
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Figure 2.4: Normalized error for fi when Su has a position error of .001" 
(solid) and .003" (dashed). The sample parameters are t T — (4.0, —j. 20) 
and \i r — (1.2,— j. 12) with d = .3". Rectangular waveguide fixture. 
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Figure 2.5: Normalized error for e when 5n has a sample length er- 

ror of +.002" (solid) and —.002" (dashed). The sample parameters are 
£ r = (4.0, — j.20) and (X r = (1.2,— . 1 2 ) with assumed d = .3". Coaxial 
airline fixture. 


13 


RE( / C6) % ERROR IM( / a) % ERROR 


8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 


oo 

1 1 I 1 i 1 1 1 J 

oo 

o 





o 







o 





o 

o 





o 

o 





o 







o — 
1 





— o 
1 

oo 

: 

i i- i 

— 

— 

oq 

o — 
1 _ 

i i i i i i i i 

— o 

_ i 


8.0 

8.5 

9.0 9.5 10.0 

FREQUENCY 

10.5 

IN 

11.0 

GHZ 

11.5 

12.0 

8.0 

8.5 

9.0 9.5 10.0 

10.5 

11.0 

11.5 

12.0 



FREQUENCY IN GHZ 


Figure 2.6: Normalized error for /z when Su has a sample length er- 

ror of +.002" (solid) and —.002" (dashed). The sample parameters are 
e T = (4.0,— J.20) and /z r = (1.2,— j. 12) with assumed d = .3". Coaxial 
airline fixture. 
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Figure 2.7: Normalized error for e when Sn has a sample length er- 

ror of +.002" (solid) and —.002" (dashed). The sample parameters are 
e P = (4.0, —j. 20) and fj. T = (1.2, —j. 12) with assumed d = .3". Rectangular 
waveguide fixture. 
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Figure 2.8: Normalized error for fi when S n has a sample length er- 
ror of +.002" (solid) and —.002" (dashed). The sample parameters are 
e r = (4.0, — j.20) and /x r = (1.2, — j.12) with assumed d = .3". Rectangular 
waveguide fixture. 
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Chapter 3 

Partially-Filled Waveguides 


Obtaining accurate e and /z values from S parameter measurements using a 
completely-filled waveguide geometry requires a tight fit between the dielec- 
tric material and the walls of the waveguide. Alternate sample geometries 
will now be explored where a known air gap is introduced between the 
material sample and one of the walls of the rectangular waveguide. Two 
different sample geometries for the dielectric sample are shown in Figure 3.1 
which are of interest. The main motivation for considering these geometries 
arises from a potentially less sensitive S parameter quantity. This sensitivity 
could introduce error in the desired material parameters and would primar- 
ily be due to air gaps between the sample and fixture walls. These air gaps 
arise from errors in sample dimensions and thermal expansion. Reducing 
the sensitivity of the S parameter measurements should facilitate the pro- 
cess of extracting accurate values for the constitutive parameters. Another 
motivation is the potential of increased flexibility in sample sizes these al- 
ternate sample geometries may provide. A greater choice of measurements 
exists because of the variety of material sample thicknesses possible. 

Compared with the completely-filled case, the numerical acquisition of 
the S parameters requires a significantly greater amount of numerical calcu- 
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lation. Because explicit relationships between the constitutive parameters 
and S parameters are not available, an iterative Newton-Raphson search is 
used to extract e and p from a set of S parameter data. The search involves 
iterating on e and n until the following two functionals are driven to zero: 

77' ncal e^ea* 

and 

F 2 = S~' - S™ a ’. (3.2) 

The calculation of the S parameters themselves involves a modal expan- 
sion of the E and H field quantities in all three regions of the waveguide and 
then matching the tangential field components at each material junction 
(z=0 and z=d ). The matching process is essentially a three dimensional 
method of moments solution using entire basis functions. 

The problem can be reduced to matching the fields at only one material 
junction (z=0) if a self-consistent technique is employed (see Appendix A). 
That is, the problem can be simplified to the semi-infinite case depicted in 
Figure 3.2 [10]. The expressions for the S parameters of a two-port system 
are then given by 

Sn = [*n] + (/ - M[P)M[P]) (■ s 2i] (3-3) 


(3.1) 


and 


s» = M[P] (i - M[P] MP])- 1 M- (3.4) 
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The variables Sij are the modal scattering matrices for the semi-infinite 
geometry and [P] is the modal propagation matrix corresponding to length 
d. 


For the one-port case, the self-consistent expression for Sn is 

Sll = [sji] -I- [s 12 pf'] {/ — [■S22][-^ / -^ , j} [$ 21 ] (3-5) 


where 

[*/'] = [ft] {M - Mlftllftl 

’ (^ + [ a n][ft][ft]) ‘M}[ftl- (3-6) 

The variables [Pi] and [P 2 ] are the modal propagation matrices correspond- 
ing to the lengths d and £, respectively. 

I Field Expansions — Geometry A 

The required elements for the scattering matrices [s*j] are obtained through 
a modal matching technique for the field expansions in both regions. The 
tangential field matching yields expansion coefficients for the expansion 
modes. The expansion modes (TE X ) for Geometry A in a rectangular fixture 
are given below* 

Region I 

E y = (e~ r ' z + a ie rii )$, + f; a n $ n e r "‘ (3.7) 

n= 2 

and 

00 

= -Y\(e- r ^ - a ie r > 2 )$ a + Y, a n Yn* n e r " z (3.8) 

n=2 
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z~0 


Figure 3.2: Waveguide regions. 






where 


*n = 

1 2 . m r 

(3.9) 

\ — sin — x, 
V a a 

r n = 

(t )’-*• Md 

(3.10) 

Y n = 

r n 

jZ 0 k 0 

(3.11) 


Region II 

Ey= £>n*ne-^ (3.12) 

n=l 

and 

oo 

= b n Y tn V n e-^ (3.13) 

n~ 1 


where 


^ n = C sin k 




*n = sin k xun {a - x ) 

7n = *«,„ - *2^, = *£,„ “ 


for 0 < ® < /, 
for t < x < a, 


Y In 

Ion 


_ Tn 

j Z 0 k 0 fs r 

— _ 2 n_ 

J ^ O^O 


for 0 < a; < f, 
for t < x < a 


(3.14) 


and C is a normalization constant. The constants k Xln and k XOn are found 
using an iterative Newton-Raphson technique and satisfy the following Iran* 
scendental equation [5] 


-k Xln cos k Xl J sin k Xl}n (a - t) = fi r k Xl)n sin k Xln t cos k XOn (a - t). (3.15) 

Continuity of the field components is now imposed at the junction z— 0. 
Next the equations are multiplied by ( q = 1, 2, TV) and integrated from 
x=0 to x—a. The final matrix form of the solution is given in Appendix B. 
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II Field Expansions — Geometry B 


For geometry B, because both TE„ and TM„ modes are excited, the expan- 
sion modes for both tangential field components, y and x, are needed to 
solve for the expansion coefficients. Below are the field expansions for this 
particular geometry. See Appendix C for the definitions of the functions 
and constants. 


Region I 


E u = 


E x = 


H x = 


Region II 


(e- r «* + a^e^AoZ^^ + £ a" e r " z A n Z^ (3.16) 

n=l 

+ a*C n Z*)*He r -, (3.17) 

n=l 

(_ e - r - + + f>*' A n + a*B n )<I>'e r »% (3.18) 


n=l 


n y = -Ewy-’ 


n— 1 


E v = 

n— 0 

H x = -b$A%*l u t-’# a 

oo 


n=l 


B, = EW 

H'i 


71—1 


(3.19) 


(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 


The constants k£f n and k^ n satisfy the following transcendental equation [5] 
and are found using an iterative Newton-Raphson technique: 


- k sin k™ t 

yin yin 


cos k* 1 ( b 

yttn V 


0 = erk“cos k*i 


M 


s ' n K!J b 


t). (3.25) 
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The constants k^ ln and k^ n satisfy the following: 


~ k L cos k viJ sin k iS b ~ 0 = ^ k L sin k mJ cos k L( b ~ 0- ( 3 - 26 ) 

As before, the continuity of the field quantities is imposed across the 
junction z= 0. The first and third equations are multiplied by 3^ ( q = 
0 , 1 ,..., AT) and the second and fourth equations by (g = 1,2, ...,7V). 
Next the equations are integrated from y— 0 to y=b. The final matrix form 
of the equations is given in Appendix C. 

Ill Calculations 

Computer software was generated to calculate Sn and S 21 for the discussed 
sample geometries, A and B. The following plots demonstrate the sensi- 
tivities for Sn and S 21 for various gap dimensions and material constants. 
Additional plots are provided to demonstrate the error in extracting e and 
\i when a gap is present in the S parameter data and processing this data 
as if no gap were present. 

The accuracy of the S parameters, of course, depends on the number of 
terms or modes kept in the field expansions. The number of modes needed 
to insure convergence is discussed in Appendix D. 

III.l Geometry A 

Figures 3.3 and 3.4 illustrate S n and S 2 i for t = .2 ",.5", .7”, and .9" in 
an X-band fixture (.4"x.9") when e r = (4.0, — j.20), /x r = (1.2, — j\12) and 
d—. 3". The variation of Sn and S 2 i decreases rapidly as the guide becomes 
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completely filled which indicates a small air gap to the short wall may not 
be that critical. 

Figures 3.5 and 3.6 illustrate Sn and S 21 for e r = 2,5,8,11, and 14 in 
an X-band fixture (.4"x.9") when fi r = (1.,— j'0.), <=.l" and d=. 3". The 
variation of Sn is quite sensitive for any loading and S 2 i would show more 
variation if loss was present in the sample. 

Figures 3.7 and 3.8 illustrate the error in extracting the material pa- 
rameters assuming the cross section of the waveguide is completely filled 
but actually having an air gap along the short wall. The air gap is .002". 
As anticipated, the error is very small or essentially zero since the electric 
field in the air gap is polarized parallel to a conducting boundary and hence 
is very small. 

Figures 3.9 and 3.10 illustrate the error in extracting the parameters 
assuming the dimension t = .490" for Geometry A but actually having an 
additional air gap of .002". The actual thickness of the material cross- 
section is t = .488". The error here is greater compared to when the 
guide is nearly completely filled. The amount of error is more significant 
in the imaginary components of the material parameters with e reaching a 
maximum of about 7%. The error in the real components remains under 
1 %. 


Note that, in general, error in the extracted parameters is increased 
where nulls occur in the S parameter data. 
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Figure 3.3: Calculated S n for t = .2", .5", .7", and .9" in an X-band fix- 
ture (A" x .9") when e T = (4.0, — j. 20), fi r = (1.2, — j. 12) and d = .3". 
Geometry A. 
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Figure 3.4: Calculated S 21 for i = .2", .5", .7", and .9" in an X-band fix- 
ture (.4" x .9") when e r = (4.0,— j. 20), \i r = (1.2, — j.12) and d, = .3". 
Geometry A. 
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Figure 3.5: Calculated Sn for e r = 2,5,8,11, and 14 in an X-band fixture 
(.4" x .9") when fi r = (1., —j 0.), t = .1" and d — .3". Geometry A. 
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Figure 3.6: Calculated £21 for £r = 2,5,8,11, and 14 in an X-band fixture 
(.4" x .9") when fi r = (1., — j’0.), t = .1" and d = .3". Geometry A. 


29 





8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 


o _ 

1 1 1 1 I 1 L 1 

1_° 

or o 
O o 
CK 
C£ 

LxJ o 
^ T 
O 

CN — 

(\) i 
^ <=> 

O 









p 

O 

O 

~~T 

p 

— CNJ 

1 

O 
— rn 
1 

p 

































— 

1 r I 1 1 T 1 HZ 1 



f n 1 1 1 1 1 1 1 r~7 

'8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12. 0 1 


FREQUENCY IN GHZ 

8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 



Figure 3.7: Normalized error in e for an .002" air gap along the short wall for 
a rectangular guide with e r = (4.0, —.7.20), fi T = (1.2, —j. 12) and d = .156". 
Geometry A but assumed completely filled. 
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Figure 3.8: Normalized error in \i for an .002" air gap along the short 
wall for a rectangular guide with e r == (4.0,— j. 20), = (1.2,— j. 12) and 

d = .156". Geometry A but assumed completely filled. 
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Figure 3.9: Normalized error in e for an additional .002" air gap along the 
short wall for a rectangular guide with e r = (4.0, —j. 20), fi r — (1.2, —j. 12), 
d = .156" and assumed t = .490". Geometry A. 
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Figure 3.10: Normalized error in fi for an additional .002" air gap along the 
short wall for a rectangular guide with e T = (4.0, — j.20), ii T = (1.2, —j. 12), 
d = .156" and assumed t — .490". Geometry A. 
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III.2 Geometry B 

In this section Geometry B will now be considered. Figures 3.11 and 3.12 
illustrate Sn and S 21 for t = .05", .1", .2", .3" and .4" in an X-band fixture 
(.4"x.9") when e r = (4.0, —j. 20), fi r = (1.2,— j.12) and d—.Z". In contrast 
to Geometry A, the variation in Sn and S 21 remains somewhat sensitive 
as the guide becomes completely filled. This indicates that a small air gap 
along the long wall may still be critical. 

Figures 3.13 and 3.14 illustrate Sn and S 2 i for e r = 2,5,8,11, and 
14 in an X-band fixture (.4"x.9") when fi T = (1.,— j'0.), <=. 1" and d~. 5". 
The variation of the S parameters for Geometry B is similar to that of 
Geometry A. Sn is very sensitive to any loading and S 2 i would vary more 
if loss were present in the sample. 

Figures 3.15 and 3.16 illustrate the error in extracting the material pa- 
rameters assuming the cross section of the waveguide is completely filled 
but actually having an air gap along the long wall. The air gap is .002". 
It is seen that the error is significantly larger than if the air gap is along 
the shorter wall. Whereas the error for an air gap along the shorter wall is 
essentially zero (see Figures 3.7 and 3.8), the air gap along the long wall for 
this case creates a maximum error of 3.5% in Im(e). This is significant in 
that when one desires to make measurements with a completely-filled guide, 
a small unintended air gap along the long wall will alter the expected S pa- 
rameter measurements which increases the difficulty in extracting accurate 
constitutive parameters. 

Figures 3.17 and 3.18 illustrate the error in extracting the parameters 
assuming t = .300" but actually having an additional air gap of .002" along 
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Figure 3.11: Calculated Su for t = .05”, .1", .2”, .3” and .4” in an X-band 
fixture (.4" x .9”) when e T = (4.0, —j. 20), fi r = (1.2, —j. 12) and d = .3". 
Geometry B. 
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Figure 3.12: Calculated £21 for < = .05", .1", .2", .3" and .4" in an X-band 
fixture (.4" x .9") when e r = (4.0, —j. 20), /i r = (1.2, — j.12) and d = .3". 
Geometry B. 
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the long wall. The actual thickness of the material was t = .298". It is seen 
that the error can be very large for both the extracted e and fi parameters, 
particularly in the imaginary components. In the lower end of the frequency 
band (under 10 GHz), the error in Re(e) is slightly less than the error when 
the guide is nearly filled while the error in fi remained about the same. It 
is anticipated that the most accurate parameter extraction for a partially- 
filled waveguide occurs when the thickness of the sample along the shorter 
wall is chosen such that the S parameters have the smallest dependency 
upon t while keeping the dimension of the sample along the longer wall 
as large as possible. In this case, the sample thickness was not chosen to 
minimize the extraction error. 
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Figure 3.13: Calculated Su for £r = 2,5,8,11, and 14 in an X-band fixture 
(.4" x .9") when fi r = (1., — j 0 . ) , t = .1" and d = .5". Geometry B. 
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Figure 3.14: Calculated 5 2 i for £*• = 2,5,8,11, and 14 in an X-band fixture 
(.4" x .9") when fi T = (1.,— jO.), t = .1" and d, = .5". Geometry B. 
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Figure 3.15: Normalized error in e for an .002" air gap along the short 
wall for a rectangular guide with e r = (4.0, — j.20), \l t = (1.2,— j. 12) and 
d = .I65' r . Geometry B but assumed completely filled. 
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Figure 3.16: Normalized error in n for an .002" air gap along the long 
wall for a rectangular guide with e r = (4.0,— j. 20), fi r = (1.2,— j. 12) and 
d — .156". Geometry B but assumed completely filled. 
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Figure 3.17: Normalized error in e for an additional .002" air gap along the 
long wall for a rectangular guide with e r = (4.0, —j. 20), fi T = (1.2, —j. 12), 
d = .156" and assumed t = .300". Geometry B. 
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Figure 3.18: Normalized error in /z for an additional .002" air gap along the 
long wall for a rectangular guide with e T = (4.0, —.7.20), fi T = (1.2, —j. 12), 
d — .156" and assumed t = .300". Geometry B. 
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Chapter 4 

Parameter Extraction 


The reflection and transmission characteristics of material samples in trans- 
mission line fixtures allow the extraction of the constitutive parameters for 
the material. As indicated in Chapter 2, simple closed-form expressions ex- 
ist when the fixture cross section is completely filled in the sample region. 
The only requirement on sample dimension is that the length should not be 
in the region of A p /2 to avoid nulls in the reflection measurement. The pres- 
ence of a null can potentially result in error due to a poor signal-to-noise 
ratio. 

The parameter extraction process is not as simple when the cross sec- 
tion of the sample region is not completely filled. An iterative procedure 
has to be initiated to acquire the constitutive parameters as indicated in 
Chapter 3. This is a result of many modes being generated to satisfy the 
required boundary conditions in the sample region. As before, the sam- 
ple dimensions should not be in the vicinity of \ g /2 to produce a null in 
the reflected field. The guide wavelength of a partially-filled fixture is also 
dependent upon the sample thickness in the cross section plane as well as 
the constitutive parameters of the sample. The sample thickness should 
be sufficient to produce Sn and S 21 parameters with comparable dynamic 
range. This occurrence will minimize the possibility of the reflected and 
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transmitted fields not sensing the material sufficiently. It would be un- 
desirable to extract constitutive parameters from measurements when the 
reflected signal is too close to zero or the transmitted signal is too close to 
unity. Reasonable choices for sample thickness and length can be obtained 
from parameter studies as illustrated in Figures 4.1 through 4.4 as well as 
Figures 3.3 through 3.6 and 3.11 through 3.14. 

Note that in Figure 4.3 the amplitude as well as the phase of the Sn 
parameter changes significantly when the waveguide is nearly completely 
filled ~ 0.4") for Geometry B. This indicates that one should avoid taking 
measurements for the completely-filled case as a small, undesired air gap 
along the longer waveguide wall would significantly alter the values of the 
S parameters. This in turn would hinder the accuracy of the extraction pro- 
cess of the constitutive parameters. As long as the air gap is large enough 
and in the range of t where both S parameters do not fluctuate greatly, 
parameter extraction ought to be accurate. Also it may be assumed that 
when taking measurements in the Geometry A configuration, significant er- 
ror in the measured S parameters may occur because of a small unintended 
air gap along the long wall. This is the case when the dielectric sample is 
not sufficiently flush against the waveguide walls. 

Geometry B may be the preferred geometry configuration provided that 
the air gap along the long wall is sufficiently large and a small error in t will 
not significantly change the S parameters. This is suggested in Figure 4.1 
where it is seen that a small air gap along the short wall does not alter the 
S parameters as drastically as in Geometry B. Then it is assumed that any 
small unintended air gap along the short wall will not significantly affect 
the S parameter measurements and accurate parameter extraction should 
occur. 
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Figure 4.1: Calculated Sn (solid) and S 21 (dashed) for varying t with 
e T = (4.0, — y.20), n r = (1.2, —j. 12), d = .3" and / = 10 GHz. Geome- 
try A. 
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Figure 4.2: Calculated Sn (solid) and S 2 i (dashed) for varying d with 
e T — (4.0, -j.20), n r = (1.2, — j.12), t = .5" and / = 10 GHz. Geome- 
try A. 
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Figure 4.3: Calculated Sn (solid) and S 2 i (dashed) for varying t with 
t T — (4.0, —j. 20), ^ = (1.2, — j. 12), d — .3" and / = 10 GHz. Geome- 
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Figure 4.4: Calculated S n (solid) and S 2 i (dashed) for varying d with 
e r = (4.0, —j. 20), fir = (1.2,-j.l2), t = .2" and / = 10 GHz. Geome- 
try B. 
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The following sections focus on the measured S parameters and ex- 
tracted constitutive parameters for two different port configurations. The 
material sample used in these examples was measured using a coaxial fix- 
ture to obtain independent values for the constitutive parameters. These 
material parameter values are shown in Figures 4.5 and 4.6. Representa- 
tive values for these parameters are e r = (4.0, —j. 20) and = (1.2, —j. 12) 
between 8 and 12 GHz. The samples were cut from a .156" thick slab of 
material and were appropriately positioned in the X-band fixture for the 
different geometry configurations. The fixture used is illustrated in Fig- 
ure 4.7. The calibration technique applied to the measured S parameters is 
discussed in [9]. Note that an additional phase factor of e~* k ’ d was applied 
to the S 21 . This accounts for the apparent increase of the fixture length 
when the sample was inserted in the fixture. This phase factor is not re- 
quired if the fixture length is physically increased by a sample holder. 


Just how accurate the extracted parameters are depends, of course, on 
the level of error terms in the S parameter measurements. Figures 4.8 
through 4.11 demonstrate how much the accuracy of the extracted consti- 
tutive parameters is affected by the level of error added to the S parameter 
data for both geometry configurations. The level of the error term added 
to the S parameters is defined to be 

ERROR LEVEL dB = 201og 10 |S 0 |ERROR LEVEL 

It is seen that the percent error for the real part of complex e and /i remains 
moderate as the error level increases. However for both geometries, the 
imaginary components of the extracted permittivity and permeability suffer 
greatly if the error level is not kept sufficiently low. Note that error levels 
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Figure 4.5: e r values using a coaxial airline fixture. 
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Figure 4.6: fi r values using a coaxial airline fixture. 
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Figure 4.7: Illustration of X-band fixture used for material measurement. 
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of —20 and —40 dB represent S parameter error percentages of 10 and 
1 percent, respectively. As seen before, a small error in the real part of 
an extracted parameter results in a large error in the imaginary part. The 
imaginary part of the material parameter is subject to larger errors because 
the real part of the material parameter is much larger than the imaginary 
part. 
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Figure 4.8: Normalized error in extracted permittivity vs. error level added 
to S parameters. Geometry A with e r — (4.0, — j.20), /u r = (1.2, — j . 1 2 ) , 
t = .490", d = .156" and / = 10 GHz in an X-band fixture. 
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Figtire 4.9: Normalized error in extracted permeability vs. error level added 
to S parameters. Geometry A with e T — (4.0, — j.20), p r — (1.2,— j. 12), 
i = .490”, d = .156” and / = 10 GHz in an X-band fixture. 
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Figure 4.10: Normalized error in extracted permittivity vs. error level added 
to S parameters. Geometry B with e r = (4.0, — j.20), n r — (1.2,— j. 12), 
t = .3", d = .156" and / = 10 GHz in an X-band fixture. 
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I Two-Port Approach 

Parameter extraction for the two-port case first requires the measurement 
of the reflected and transmitted field quantities (i.e., Sn and S 21 ) for a 
sample-filled waveguide. After a calibration has been performed on the 
measured data, the constitutive parameters for the sample may now be 
obtained via an iterative two-dimensional Newton-Raphson technique. 

Figures 4.12 and 4.13 show the measured Sn and S 21 data for Geome- 
try A with sample dimensions t = .490" and d = .156". For comparison the 
S parameters, calculated using the constitutive parameters of the coaxial 
line fixture, are also plotted. It is seen that the calculated and measured 
values are very similar for both S parameters except after about 11 GHz. 
This indicates that the parameter extraction process should produce accu- 
rate e and /r values at frequencies below 11 GHz. 

Figures 4.14 and 4.15 illustrate the measured and calculated S parame- 
ter data for Geometry B with sample dimensions t = .300" and d — .156". 
Notice that around a frequency of 10.5 GHz the measured Sn begins to 
diverge considerably from the expected calculated Sn. Also after 11 GHz 
the measured S 21 begins to differ significantly from the calculated values of 
S 21 • These discrepancies may be explained by the fact that there was not 
sufficient contact between the sample material and the wall of the waveg- 
uide. Also error in the calibration technique may have occurred. In any 
case parameter extraction from these S parameter measurements for Ge- 
ometry B cannot be expected to be accurate after around 10.5 GHz for this 
particular case. 

This Newton-Raphson method uses two functionals, F\ and F 2 , which 
consist of calculated and measured values of S n and S 2 i and are defined as 
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F t = s;f ( C ,tt) - sz’°‘ 


(4.1) 


and 


F, = SJf {',?) - SJ~. 


(4.2) 


An initial set of values for and m are used to find the calculated Sn 
and S 21 parameters. The next set of e,-+i and m+i are found using 


e»+i 


f*i + 1 

where 


A ti 


Am 


and 


J - 


Ci “I - AC{ j 

m + 





J 





dF\ BFx 

dc dfj 


dF* dFi 
6 c B11 


(4.3) 

(4.4) 


(4.5) 

(4.6) 


(4.7) 


The iteration process is continued until both functionals are driven to 
zero. The extracted set of constitutive parameters is then the values used 
in the final iteration. 

Figures 4.16 and 4.17 show the extracted e r and m values, respec- 
tively, for Geometry A from the set of measured data given in Figures 4.12 
and 4.13. As expected because of the reasonable measured S parameter 
values, the extracted values for e and /x are fairly stable and accurate below 
11 GHz, especially the real components. The imaginary components are 
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not as accurate as the real parts, but it was shown earlier that the accuracy 
of the imaginary components would be affected more significantly by noise. 


Figures 4.18 and 4.19 show the extracted e r and /x r values for Geome- 
try B from the set of measured data given in Figures 4.14 and 4.15. Al- 
though the values for extracted e r are a little low, the values for extracted 
\i T are accurate throughout the frequency spectrum. This is true despite 
the S parameters diverging from the calculated values after 11 GHz. 

Parameter extraction seems to have worked fairly well for both geome- 
tries. The real part of the extracted values were around what they were 
expected to be. But there was some inaccuracy in the imaginary parts of 
the extracted constitutive parameters. This inaccuracy is a result of the 
real part of the extracted values being much larger than the imaginary part. 
Vectorially, a small angular change of the complex extracted value results 
in essentially the same value for the real part (|e|cosa ~ |e| for small a) 
while the imaginary part can experience dramatic changes (|e| sin a % |er|o). 
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Figure 4.12: Measured (solid) and calculated (dashed) Sn for t = .490" and 
d = .156" in an X-band fixture (.4"x.9"). Calculations use e r = (4.0, —j. 20) 
and fi r = (1.2, —j. 12). Geometry A. 
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Figure 4.13: Measured (solid) and calculated (dashed) S 21 for t = .490" and 
d = .156" in an X-band fixture (.4"x.9"). Calculations use e r = (4.0, —.7.20) 
and /i r = (1.2,— j. 12). Geometry A. 
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Figure 4.14: Measured (solid) and calculated (dashed) Sn for t = .300" and 
d = .150" in an X-band fixture (.4"x.9"). Calculations use e r = (4.0,— y.20) 
and fi r — (1.2,— j. 12). Geometry B. 
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Figure 4.15: Measured (solid) and calculated (dashed) S 21 for t = .300" and 
d = .156" in an X-band fixture (.4"x.9"). Calculations use e T — (4.0, —j. 20) 
and fi T = (1.2,— j. 12). Geometry B. 
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Figure 4.16: Extracted e r for Geometry A with t = .490" and d .156" in 
an X-band fixture (.4"x.9"). 
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Figure 4.17: Extracted // r for Geometry A with t = .490" and d = .156" in 
an X-band fixture (.4"x.9"). 
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Figure 4.18: Extracted e T for Geometry B with t = .300" and d — .156" in 
an X-band fixture (.4"x.9"). 
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Figure 4.19: Extracted /x r for Geometry B with t = .300" and d = .156" in 
an X-band fixture (.4"x.9"). 
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Figure 4.20: Single-port fixture geometry. 

II Single-Port Approach 

The necessary S parameters for material parameter extraction can be ac- 
quired from a two-port fixture. However, a single-port fixture may be used 
to simplify the measurement process and requirements of fixture design at 
elevated temperatures. For this simplification, additional signal processing 
is required to recover the Sn and S 2 i parameters. Figure 4.20 illustrates 
the single-port fixture geometry. The reflected field for this geometry can 
be expressed as 


E T = S„ 


c2 -j20„l 
J 21 e 


1 + Sne-'W 


(4.8) 


when only the dominant mode can propagate in the air-filled portion of the 
guide and /3 0 is the propagation constant for that mode. It is important 
to account for loss in the propagation constant, /3 0 , if the fixture itself has 
sufficient loss. The self-consistent expression for this fixture geometry is 
given in Appendix A which accounts for all modes in the air-filled portion 
of the fixture. The other higher order modes are not of great importance 
when the sample/short separation distance, £, is large. 


Equation (4.8) contains the S parameter information required to de- 
termine the material parameters. This expression is simply the geometric 
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summation of all the multiple reflected and transmitted signals due to the 
short in the fixture. Figure 4.21 illustrates the terms of this summation 
in a time domain plot for a completely-filled waveguide cross section and 
d — .156". The first signal (f=0) in this figure is simply S n while the 
second is — These desired quantities can be recovered from the total 
reflected field through signal processing techniques. One technique is time 
domain gating [11] which may be undesirable due to the small single mode 
bandwidth available in waveguide fixtures. An alternative is to use a mov- 
ing average technique upon the frequency domain data. This technique is 
conceptually comparable to the time domain gating approach because it du- 
plicates a discrete convolution. This particular moving average performs a 
weighted average around the t 4 * 1 point using the neighboring n points. The 
greater n becomes, the tighter the effective time gating window becomes. 

The Sn parameter can be recovered from the measured data by simply 
applying a moving average to the data when the phase reference is placed at 
i — 0 and using the proper number of neighboring points. The S 2 i parameter 
is recovered in a similar manner but with the data initially phase adjusted 
to place the second signal at t= 0. This is accomplished by multiplying the 
measured data by the phase factor of 

The performance of this technique is demonstrated in Figures 4.22 
through 4.27 where the distance from the front of the sample to the short 
( d + i) was 10.012" and an X-band fixture was used. Figures 4.22 and 4.23 
contain the recovered S parameters for a completely-filled cross-sectional 
waveguide sample region. Figures 4.24 through 4.27 contain the recovered 
S parameters for the partially-filled case for Geometries A and B. Because 
of the type of processing that was used, the recovery of the S parame- 
ters at the frequency band edges was not good. Note that for the moving 
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average approach to succeed in recovering the desired S parameters, the 
sample/short separation distance, £, has to be sufficiently large. A length 
of at least six inches is required. Also note that the sample length should 
be less than A ff /2 to avoid nulls in the S parameters which is a similar re- 
quirement when the cross-section is completely filled. This requirement is 
based not only on the need to maintain a good signal to noise ratio but 
also to avoid nulls since the moving average technique will “fill” them in. 

Now that the Sn and S 2 i data has been determined, the extraction 
procedure of the two-port method is used to find the constitutive parame- 
ters. Figures 4.28 through 4.33 show the extracted e and fi values for the 
completely-filled case and the two geometries of Figures 4.22 through 4.27. 
Note that the measured and calculated values of the S parameters are in 
fairly good agreement. This indicates that accurate constitutive parameter 
extraction should occur. As expected, Figures 4.28 through 4.33 suggest 
that accurate parameter extraction has indeed been achieved at least for the 
real components. The imaginary components suffered a bit of inaccuracy 
due mostly to the small amount of error introduced during the process of re- 
covering the S parameters from the one-port data. For the completely-filled 
case and both geometry configurations, it was found that because of this 
error the imaginary component of e tended to be low while the imaginary 
part of fi tended to be high. 
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Figure 4.22: Recovered Sn from a measured one-port file (solid) com- 
pared with Sn from a two-port measurement (dashed). Completely-filled 
cross-section with d = .156" in an X-band fixture. 
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Figure 4.23: Recovered S 21 from a measured one-port file (solid) com- 
pared with S 21 from a two-port measurement (dashed). Completely-filled 
cross-section with d = .156" in an X-band fixture. 
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Figure 4.24: Recovered Sn from a measured one-port file (solid) com- 
pared with Sn from a two-port measurement (dashed). Geometry A with 
i = .490" and d = .156" in an X-band fixture. 
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Figure 4.25: Recovered S 21 from a measured one-port file (solid) com- 
pared with S 21 from a two-port measurement (dashed). Geometry A with 
t — .490" and d = .156" in an X-band fixture. 
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Figure 4.26: Recovered Sn from a measured one-port file (solid) com- 
pared with Si i from a two-port measurement (dashed). Geometry B with 
t = .300" and d — .156" in an X-band fixture. 
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Figure 4.27: Recovered S 21 from a measured one-port file (solid) com- 
pared with S 2 i from a two-port measurement (dashed). Geometry B with 
t = .300" and d = .156" in an X-band fixture. 
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Figure 4.28: Extracted e r from a one-port file for a completely-filled waveg 
uide with d = .156" in an X-band fixture. 
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Figure 4.29: Extracted fx r from a one-port file for a completely-filled waveg- 
uide with d = .156" in an X-band fixture. 
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Figure 4.32: Extracted e T from a one-port file. Geometry B with i 
and d = .156" in an X-band fixture. 
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Figure 4.33: Extracted /z r from a one-port file. Geometry B with t 
and d — .156" in an X-band fixture. 


m 



Chapter 5 
Conclusions 


Waveguide structures can be used to ascertain the constitutive parameters 
e and fi of a material. Traditionally for a rectangular waveguide fixture its 
cross section is completely filled with the material to be measured. But with 
a high temperature environment, air gaps can appear which cause error in 
the S parameter measurements. This lessens the accuracy of the constitu- 
tive parameters that are extracted from these S parameter measurements. 
This report examined the effect of air gaps on the accuracy of the extracted 
parameters and provides alternate approaches to account for these air gaps 
in the extraction process. An added benefit of introducing an air gap be- 
tween the fixture’s side walls and sample allows a greater choice in sample 
dimensions to provide for a desired level of measured signals. Also included 
is a description for a single-port approach to obtain the necessary quanti- 
ties (Sn and S21) obtained from the more traditional two-port approach. A 
single-port approach has mechanical benefits for high temperature heating 
considerations. 

Errors in the measurement originate from several sources. This report 
focused upon the effect of air gaps but also considered sample positioning 
and dimensions error as well. Errors from calibration or fixture connections 
were not considered. 
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The effect of air gaps along both side walls of the fixture was exam- 
ined. A small air gap oriented along the short wall of a rectangular waveg- 
uide (Geometry A) was shown to have little effect on the S parameters 
which means that parameter extraction should remain accurate. In con- 
trast, a small air gap oriented along the longer wall of a rectangular waveg- 
uide (Geometry B) alters the S parameters significantly so that the recovery 
of an accurate set of constitutive parameters would become difficult. But if 
the size of the air gap in Geometry B is made larger and known beforehand 
the S parameters may not change very much with small changes in the size 
of the air gap. Thus the preferred orientation of the sample in the rect- 
angular waveguide is Geometry B because small gaps along the short wall 
that may appear with high temperature variation does not seem to impair 
accurate parameter extraction as drastically as small gaps along the long 
wall. 

Other examined errors which occur in any rectangular waveguide mea- 
surements are those due to phase offset error in Sn and sample length error. 
The phase offset error which is caused by the sample not being placed where 
the calibration reference was placed effected a large amount of error in the 
imaginary part of the extracted constitutive parameters. Thus the posi- 
tioning of the material sample must be fairly accurate. On the other hand 
the accuracy of the extracted parameters is not nearly as sensitive to the 
error between the actual and assumed lengths of the material. 

The two-port and single-port approaches to S parameter measurement 
were discussed. The two-port approach requires simply acquiring Sn and S 2 i 
and initiating the extraction process. The one-port approach would be more 
suitable to high temperature variation since only one end of the waveguide 
would be necessary to acquire Sn. The one-port method involves recover- 
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ing the equivalent Sn and S21 parameters of a two-port system from the 
one-port Sn measurements. Because the one-port approach requires an ad- 
ditional signal processing step, it does not yield as accurate e and // values 
as in the two-port approach. 

The accuracy of the presented single-port technique relies heavily upon 
knowing the separation distance between the material sample and the short 
very accurately. This may have limitations for elevated temperature mea- 
surements. An alternative to this technique but still using a single-port fix- 
ture is to measure the return from two material samples of different lengths 
which have their far side positioned against a short. This technique is simi- 
lar to the presented technique but with no separation distance between the 
material sample and the short. With this alternate technique, two different 
sample lengths have to be measured to acquire two independent responses. 
The presented technique acquires two independent responses through signal 
processing. Such processing is not possible when the separation distance is 
zero. 

Finally, significant error is introduced in the imaginary component of 
the extracted e and n values when the loss is small. This result occurs in 
any transmission line measurement. Vectorially examining the extracted 
values indicates the explanation. A small change in the magnitude or angle 
of this vector will result in only a small change in its real component while 
a significant change will occur in its imaginary when the loss is small. 
However, the uncertainty in the loss value is minimized when a sufficient 
sample length is used. 



Appendix A 


The Self-Consistent Technique 


The self-consistent procedure is used to simplify the analysis of the two- 
and single-port waveguide problems. Potentially, one method used to find 
the reflected and transmitted fields in a two-port problem is to expand the 
fields in the three different regions of the waveguide and then enforce the 
boundary conditions at the two material interfaces. The desired solution 
then requires a matrix inversion. The size of the required matrix to be 
inverted can be reduced to one half the size if only one interface boundary 
(semi-infinite problem) needs to be enforced. The self-consistent technique 
provides this smaller dimensioned matrix to be inverted. This benefit re- 
sults in a faster solution when the order of the matrix becomes large. The 
desired two-interface problem can be solved with the one interface results. 
Formally the size of these matrices is infinitely large, but in practice there 
are N x N matrices where N is the number of terms kept in the modal 
expansions. 

Once the scattering matrices are known, they can be used in the formu- 
las for the reflected and transmitted fields. These formulas will be derived 
next using the self-consistent procedure for the two-port and single-port 
formulations. 
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I Two-Port Formulation 


Referring to Figure A.l, let [A + ] denote the vector containing the ampli- 
tudes of the modes of Region II traveling in the +n (- z ) direction and [A - ] 
the amplitudes of the modes traveling in the -ra ( +z ) direction. Then the 
reflected and transmitted fields can be decomposed as 

[£'] = (»»][£*] + M[A + ] (A.i) 

and 


l&) = Mp}[a-} 


(A.2) 


where [£’*] is the incident field vector and [P] is the diagonal matrix account- 
ing for the propagation phase delay of the modes over the path length d. 



[A+] 

[A-] 


Figure A.l: Waveguide regions for the two-port approach. 
The vectors [A f ] and [A - ] can be expressed as 


[A + l = [P][.„l[Pp-] 


(A.3) 
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and 


[A ] — [<S2i ] [-£'*] + [^ 22 ] [ ] • (A.4) 

Inserting Eq. (A.3) in (A.4) gives 

[A ] = [s 21 ][E‘] + [s 22 ][P][s 22 ][P][A - ]. (A-5) 

Rearranging Eq. (A. 5) leads to 

[AT] = (/ - [•aJ[P][*»][P])- 1 (A. 6) 

Substituting Eq. (A. 6) back into (A.3) gives 

[A + ] = [P]M[P] (I - MP1MP])- 1 (A -7) 


Inserting Eq. (A. 7) into (A.l) gives the equation for the reflected field and 
is given by 

[£'] = {[«„! + [.JflHIPl (/ - MPIMIPir’ [»,,]} [£’] (A.s) 

and inserting Eq. (A. 6) into Eq. (A. 2), the expression for the transmitted 
field is 

[£'] = [.„][P| (/ - MP1MP])- 1 [.„][£•']. (A.9) 

II Single-Port Formulation 

For the single-port method where a short is placed a distance £ from the 
end of the dielectric material, the self-consistent procedure is used again 
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to determine the reflected field from the scattering matrices. Referring to 
Figure A. 2, let the vectors [A* ] and [Ajj"] denote the amplitudes of the modes 
of Region II at z — 0 traveling in the +n and -n directions, respectively, and 
let [Aj] and [A^] represent the amplitudes of the modes of Region III at 
z=d traveling in the +n and -n directions. Then the reflected field can be 
expressed as 


m = [»..][£'] + [»1>)K) 

(A.10) 

where 

W) = [fllWie.Pr] + [A1K1IAJ-], 

(All) 

Mil = [ a 2l]M‘] + [ 5 22 ]Mi + ], 

(A.12) 

Ma] = [^JM^JMa] and 

(A. 13) 

M2] = M(Pi]Mr) + MM2] 

(A.14) 

where [F*] is the incident field vector, [ar] = — I is 

the scattering matrix for 


the shorted end, [Pj] is the diagonal matrix accounting for the propagation 
phase delay of the modes over the path length d and [P 2 ] is the diagonal 
propagation matrix over the path length L 


Inserting Eq. (A. 13) into Eq. (A. 14) and rearranging gives 

[A,-] = (/ - [»i.I(e,][» r ](e,l )" 1 MP.pr)- (A.is) 

Substituting Eq. (A. 15) back into (A. 13) 


\Ai\ = [^JM[^ 2 ] (J- - M(P,Pr]. (A. 16) 


Putting Eq. (A. 16) in (A. 11) 

\At ) = {[niwini 
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SHORT 


\ 



Figure A. 2: Waveguide regions for the one-port approach. 

+ [j>,i[j„)[iyi.r][jy (/ - [-.nHfiiHift])- 1 i.»][i’.]} Mrl 

= \M\\A-,) (A. 17) 

where 

[A/] = [Pil {[»»] + WlftlWIft] 

• (/ - MIM'rMr 1 M} !«]■ (a. 18) 

Inserting Eq. (A. 17) into (A. 12) and rearranging leads to 

Mr] = P - WW)" Ml^l- (A-i9) 

Substituting Eq. (A. 19) back into Eq. (A. 17) gives 

[At] = [ M ] {/ - MA/]}" 1 [i 2 ,][F]. (A. 20) 

Finally using Eq. (A. 20) in (A. 10) and the fact that [sr]=~^> the reflected 
field is 


93 




[Pi = {[»„] + [ , U ][M'] { I - MIA /']}' 1 M} |£‘). 

where 

\ m '] = M {[.»]- [.,.][ft]|/y 


(A.21) 


(A. 22) 
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Appendix B 

Scattering Matrix Formulation 
for Geometry A 


This section deals with attaining the elements of the scattering matrices 
which are to be used in the self-consistent procedure shown in the previous 
appendix for Geometry A. 


First the scattered fields are defined below for Regions I (z < 0) and 
II (z > 0) where the scattered fields denoted by a superscript s are prop- 
agating or attenuating in the - z direction for Region I and propagating or 
attenuating in the +z direction for Region II. 


Region I 



X! an$r,e r 

n= 1 


and 


= £ a n y n * n e r - 

71 rr 1 

where 



(B.l) 


(B.2) 


(B.3) 
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(B.4) 

(B.5) 


T 2 n 

Yn 


(t) ~ kl and 

r n 

jZ 0 k 0 ' 


Region II 

OO 

^’"= ne ' 7 " 1 

n=l 

and 


H '/ 1 = -E b nYin* n« _V 

n = l 

where 


(B.6) 


(B.7) 


= Csinfc xin x 
Wn ” ^ »>n fct nn (a — t) 


sinfc Xon (a - x) 


7n = kl ln - kle rf i r 


- U 2 


— fc 2 

K oi 


Yi n 

y 0n 


_ 2a 

jZ 0 koH T 

— 7n 
jZ 0 k 0 


for 0 < x < t, 
for t < x < a, 

for 0 < x < t, 
for t < x < a 


(B.8) 


and C is a normalization constant. The constants fc Xln and fc XOn satisfy the 
following transcendental equation and are found using an iterative Newton- 
Raphson technique. 


-k Xin cosk Xln t sin k x<>n (a - t) = fi r k XOn sin k Xln t cosk K0n (a-t) (B.9) 

To get the first column of the s u and s 21 scattering matrices, it is as- 
sumed that the TEi 0 mode in Region I is the excitation mode propagating 
in the +z direction and incident on the material junction (z=0). The ex- 
citation field quantities of Region I denoted by the superscript i are given 

by 
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(B.10) 


El = t,e 


-fiz 


and 


HI 


Y&e 


-Viz 


(B.ll) 


Matching the total tangential E and H fields (incident+scattered) at 
z = 0 gives 

(1 + <Zi)$i + ^ = (B.12) 

n= 2 m = l 

and 


oo oo 

-(1 - ai )Ei*i + £ a n F n * n = - Y, (B.13) 

n—2 m = l 

Multiplying the above equations by the test functions $ q (g — 1,2,3, ...), 
integrating from x=0 to x=a, keeping the first N terms of the series and 
rearranging gives 


N 


^ ' ( ”1“ ^n-^qn) 


n=l 


and 


N 


Yi6 q i = ^(a n y n <5 qn + 6 n y 0n (5gn) 


n— 1 


where 


Sqn = < 


qn 


1 for q = n 
0 for q ^ n 


= f 

Jo 


Qq'tyndz and 


Qqn = — f $qVn<Lx + / ^ q ^ n dx. 
Hr JO Jt 
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(B.14) 


(B.15) 


(B.16) 

(B.17) 

(B.18) 


In matrix form Eqs. (B.14) and (B.15) are 


[X][S] = [E] 


(B.19) 


where the matrices X,S,and E are illustrated in Figure B.l. 


1 0 0 • • • 0 

0 -1 0 • • • 

* * 

• • 

• m 

0 • * • 0-1 

Pu Pl2 • • • PiH 

P 21 p 22 . . . p 2N 

■ • • 

* 

• • 

P N1 P N2 * ’ ’ ’PNN 

Y x 0 0 ... 0 

^oAl' 1 T)2Ql2* * * XjN^lN 

0 y 2 0 • • • 0 

• * 

• f 

Y 01 Q ll Y 02Q22* * * ^OiAn 

9 9 

9 ^ 

• 

0 • • • ’0 ^ 

^>1®N1^>2^I2' ' * Yqh^JN 


[S] = 


1st COLUMN 
OF S u MATRIX 


hi 


[E] = 


. 1st COLUMN 
>OF s 21 MATRIX 



0 

Yi 

0 



1st COLUMN 
> OF EXCITATION 
MATRIX 




Figure B.l: Matrix solution for the first column of the sn and s 2 i scattering 
matrices. Geometry A. 


To get the k th ( k — 2,3 , ...,7V) columns of the Sn and s 21 scattering 
matrices, the incident fields of Region I are 


E'v 


^ e - r * 2 


(B.20) 


and 
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(B.21) 


Hi = 

Then the same procedure shown above is used. 

To get the k th columns of the S 22 and S 12 scattering matrices, the incident 
field now propagates (or attenuates) in the -z direction in Region II and 
one uses the following incident field expressions for Region II: 

£* = tf fc e^ 2 (B.22) 

and 

Hi =Y ik * k e'" z . (B.23) 

The entire matrix structure is illustrated in Figure B.2 where the X 
matrix in Figure B.2 has been partitioned into four parts. 


[S u ]| [S i2 ] 

- —1 


CNJ 



r-i 
1 — 1 

1 

“I 

[-XjJ | [-xj 

[S21] [s 22 ]_ 


[Xailj t x 22 ] 


[X 21 ]i [x 22 ] 


Figure B.2: Complete matrix solution for all four scattering matrices. Ge- 
ometry A. 
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Appendix C 

Scattering Matrix Formulation 
for Geometry B 


This appendix shows how the elements of the scattering matrices needed 
for the self-consistent procedure are attained for Geometry B. 


First the scattered fields are defined below for Regions I and II where the 
scattered fields denoted by a superscript s are propagating or attenuating 
in the - z direction for Region I and propagating or attenuating in the +z 
direction for Region II. 

Region I 



II 

* 5* 

00 

n=rO 

(C.l) 

E’ 1 = 

f + of C.zf )#”e r --, 

n= 1 

(C.2) 

H' 1 = 

a*%$'e r " + f>"4„ + 

n~ 1 

(C.3) 

** 

to 

II 

- t a nCn^e^ 

Tl— 1 

(C.4) 

where 



& 
a *■* 

II 

— ^ 

f \/T sin ^ for n = 0 

l y/fsin^cos 1 ^ for n = 1,2, ... 

(C.5) 
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* n 

r! 


z m 


z* 


Ar 

B T 


2 7rx . mry 

V6 cos T s,n_ T’ 


weoI\, 

vHoT r 


> [*s - (f ) J ' 


jwoT n , 
a 

and 


n 7T 


afe 


7T£ 


= sin— 


a 


7 TX 


cos fcj ln y 

cosfel, t 


for 0 < y < t 


co. fc t^ COsfc U 6 -J/) for<<y<6 


¥? t = cos — 


sin fcj ln y 

sinkl, f 


for 0 < y < < 


l fori<y<b 


\ 2 

r n - 


k L 2 - k i + (J) . 


(C.6) 

(C.7) 

(C.8) 

(C.9) 

(C.10) 

(C.ll) 


c n = 


(C.12) 

Region II 



E v - 

n— 0 

(C-13) 

E‘f = 

oo 

nr- 1 

(0.14) 

H' x u = 



- EWAfVire-#' + l£*«e-*«), 

n=l 

(C.15) 

H’ y n = 

y^cE^nEe-y^ 

✓ n n n 

n=l P* 

(0.16) 

where 



i = 

f 1 for 0 < y < f 
| 0 for t < y < b 

(0.17) 


(C.18) 

(C.19) 

(C.20) 
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(C.21) 



(!) : 


fn 


U>£;7, 


M 


* [(?)■- 

• M 

JWlln » 


Tif 2 


yin 


and 


a 



(C.22) 

(C.23) 

(C.24) 

(C.25) 


The constants kjjf n and kj}* n satisfy the following transcendental equation 
and are found using an iterative Newton-Raphson technique. 


sin k i!J cos KlS b - 0 = e *KL cos k mJ sin K f J b - *)■ (c.26) 

The constants k yin and fc^ )n satisfy the following: 

- k L cosk ZJ sink L( b - = ^ k L sin k L f cosfc L( 6 -0- (c.27) 

To determine the first column elements of the and S 21 matrices, 
the TE 10 mode (or TM yi0 ) in Region I is assumed to be the excitation 
mode propagating in the +z direction and incident on the material junc- 
tion (z=0). The excitation quantities exist only in Region I and are 


E' y = 

(C.28) 

Ei = 0, 

(C.29) 

Hi = -A 0 ^ 0 e- Uz 

(C.30) 
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and 



(0.31) 


Matching the total tangential E and H fields (incident -f-scattered) at 
2=0 gives 


(1 +<tf)A 0 Zi’4l+'E*“A„Z?i' n = E^ A «Z^¥”, (C.32) 

n= 1 m=0 

+ 4 C„ZZ)iH = f 

m=l 

+‘£<£25'*“) , 

(C.33) 


n— 1 


(1 - + Z«A n + ^B„)K = 

n=l 

-C-OJ" - E(tJMi?*i" + £bZK e ) 


m — 1 


and 


-E 




n=l 


lEriE^ 1 Hi HE 
m = I ft* 


= E 


(C.34) 


(C.35) 


Multiplying Eqs. (C.32) and (C.34) by the test functions 
(? = 0,1,.. .,JV), multiplying Eqs. (C.33) and (C.35) by the test func- 
tions (g ; = 1,2,..., JV ) , integrating from y=0 to y—b, keeping the first 
N terms of the series and some rearranging gives 


A 0 Z 0 6 q o 

0 


E(-“n A„Z?6, n + b" AfZ ? 1 (?"'), 

n=0 

E {-<«?*. + af C„Z*)Vn + 

1 

+6fC*Zf F,7,f } , 


n~ 1 


iTJIf 


(C.36) 


(C.37) 
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Ao6 g o = 

+ b“A*'p;” 

oo 

+ E + *ZB„)6„ + 6 + 4fB* J*f } 

n=l 

and 


0 = E(«? C "Vn 


n— 1 


+ £<?<$?) 


where 


^gn 

p/ f 

gn 

pilf 

gn 

<?;r 

Q™ 


1 for q — n 
< ? 

0 for q ^ n 

jf *f *f«e«y, 

jy^dy + erj'z^dy and 

jf + ^ f t * n q y« E d y . 


(C.38) 

(C.39) 

(C.40) 

(C.41) 

(C.42) 

(C.43) 

(C.44) 


Eqs. (C. 36), (C. 37), (C.38) and (C.39) are given in matrix form as 


[X][SJ = [E] 


(C.45) 


where the matrices X, S, and E are illustrated in Figure C.l. 

To get the k th ( k = 1,2, ...,7V) columns of the s n and s 2 i scattering 
matrices, the incident fields used in Region I are 


E‘ = A t z!'H*- r, ‘, (C.46) 

K = B k i"e~ rk ‘, (C.47) 
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(C.48) 


Hi = -A k *[e- T " 


and 


HI = 0. (C.49) 


Then the same procedure shown above is used. 


To get the k' th (k 1 = k + N) columns of the the s n 

and s 2 i scattering 

matrices, the incident fields in Region I are 


Ei = 0, 

(C.50) 

Ei = 

(C.51) 

Hi = -B k .$[,e- r »' z 

(C.52) 

and 


W y = C k >*Z,e- r *' z . 

(C.53) 

To get the k th (k = 0 , 1 ,..., AT) columns of the s 2 2 

and sn scattering 

matrices, the incident field now propagates (or attenuates) in the -z direc- 
tion in Region II, and one uses the following incident field expressions for 

Region II: 


Ei = 

(C.54) 

f for k ^ 0 

E l x = 

0 for k — 0 

(C.55) 
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(C.56) 


Hi = A¥*l M e* 


and 


H' v = 0. 

To get the k ,th (k 1 = k -f N,k = 1,2,..., N) columns of the the 
S 21 scattering matrices, the incident fields in Region I are 

K = o» 

El = 

and 

H' x = B*'Zl?e y * z 


and 


Hi = 


The entire matrix structure is illustrated in Figure C.2. 


(C.57) 
sn and 

(C.58) 

(C.59) 

(C.60) 

(C.61) 
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Figure C.2: Complete matrix solution for all four scattering matrices. Ge- 
ometry B. 
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Appendix D 

Determination of Gamma 
Constants 


The calculation of the S parameters for a partially-filled waveguide requires 
finding a set of gamma (propagation) constants. For Geometry A the 
gamma constants *y n (n = 1 , 2 , 3 ,. . .) used in the series expansions for the 
S parameters satisfy the transcendental equation 

cos k„J sin k,„(a - t) + sin cos - () = 0 

(D.l) 

where 


7n — k h n k o e rf l r — k 




(D.2) 


For Geometry B, the TM y gamma constants 7 * f (m = 0,1,2,...) must 
satisfy 

sin k"j cos k"jb - <) + cos k "j sin *=“ J b ~ l ) = 0 

(D.3) 

and the TE y gammas 7 ® (n = 1,2,3,...) satisfy 

k L cos k viJ sin k L( b “ 0 + #*■■*£„ «n k^J cos fc£ n (& - <) = 0 

(D.4) 
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where 


= *i.’ -*?+(!)’ • < D - 5 ) 

i = ( 1,OT °^‘. (D. 6 ) 

0 for t < y < b 

Each of the above transcendental equations contains an infinite number 
of roots. The method of obtaining the gamma constant for a particular 
index number n involves first finding the value which satisfies the above 
transcendental equations for the completely-filled (or empty) case. The 
exact gamma value of a completely-filled guide for Geometry A and 
index n is given by 



jk 0s /e r n r 


1 


1 - 


_A|_ 

4e r /z r 



(D.7) 


The exact gamma value for Geometry B of a completely-filled guide is given 
by 



jk 0 ^e r fi T 



(D. 8 ) 


Then the dimension t in Eqs. (D.l) through (D.4) is decreased from 
the completely-filled case by an increment At, and the initial gamma 7 ^ 
calculated in Eq. (D.7) or Eq.(D. 8 ) is used in the first iteration to find 
the gamma which satisfies the appropriate transcendental equation. This 
gamma value is then the initial value used in the next iteration where the 
current value of t in the transcendental equation is again decreased by At. 
This process is continued in increments of At until the desired dimension t 
is reached. 
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For each t used in the above procedure a Newton- Raphson iteration pro- 
cess is used to find the gamrna constants which satisfy their transcendental 
equations. 

Problems may occur in the search when At is not small enough. Also, 
finding a wrong gamma may occur when a gamma passes through zero in 
the complex plane during the iteration. A plot of how the gammas for 
different modal indices vary with t is given in Figure D.l for Geometry A 
and Figures D.2 and D.3 for Geometry B. A plot for varying frequency is 
also given in Figure D.4 for Geometry A only. 

The number of gamma constants actually needed to calculate the S pa- 
rameters depends on the desired accuracy of the S parameters. Figures D.5 
and D.6 illustrate how the S parameters vary with the number of gamma 
constants (modes) used in the calculations for Geometry A. Similarly re- 
sults for Geometry B are presented in Figures D.7 and D.8 where a number 
N modes used refers to (N + 1) TM y modes and N TE y modes. The above 
figures used e T = (4.0, —j-2) and fi r = (1.2, —j. 12) in an X-band fixture. 
The typical number of modes used to generate the plots throughout this 
report was six for both geometries. 
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Figure D.l: Plot of alpha vs. beta from t = 0 to t = a where 7 = a + jfl- 
Geometry A in an X-band fixture (.4"x.9") with e T = (4.0, -j.20), 
H r = (1.2, — j. 12) , d = 0.3" and / = 10 GHz. 


112 






Figure D.2: Plot of alpha vs. 
where 7 = 0 + j(3. Geometi 
e T = (4.0, — j.20), n r = (1.2, -j. 




Figure D.3: Plot of alpha vs. beta for TE V modes from t = 0 to t = b 
where 7 = a + j(3. Geometry B in an X-band fixture (A"x. 9") with 
e r = (4.0, -j.20), fi r = (1.2, — j.12), d = 0.3" and / = 10 GHz. 
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Figure D.4: Plot of alpha vs. beta from / = 8 GHz to / = 12 GHz 
where "y = a -f j/?. Geometry A in an X-band fixture (.4 z . 9 ) with 
e r = (4.0, — j.20), fi r = (1.2,— j.12), d= 0.3" and t = 0.6". 
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# MODES USED 


Figure D.5: Sn vs. number of modes used in calculations. Geometry A 
with t = .514", d = .3" and / = 10 GHz. 
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- 3.90 - 3.70 - 3.50 - 180 . - 90 . 






MAGNITUDE IN DB PHASE IN DEGREES 



# MODES USED 


m 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 



# MODES USED 


Figure D.6: S 21 vs. number of modes used in calculations. Geometry A 
with / = .514", d = .3" and / = 10 GHz. 
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- 5.70 - 5.65 - 5.60 - 5.55 - 5.50 - 5.45 - 180 . - 90 . 


MAGNITUDE IN DB PHASE IN DEGREES 

27 -5.24 -5.21 -5.18 -5.15 -180. -90. 0. 90. 



# MODES USED 


Figure D.7: Sn vs. number of modes used in calculations. Geometry B 
with t = .3", d = .3" and / = 10 GHz. 
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